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Abstract. We prove the existence of quasi-stationary symmetric solutions 
with exactly n > zeros and uniqueness for n = for the Schrodinger-Newton 
model in one dimension and in two dimensions along with an angular momen- 
tum m > 0. Our result is based on an analysis of the corresponding system of 
second-order differential equations. 



1. Introduction 
We consider the Schrodinger-Newton equations 

(1.1) itp t + A V> -7$V = 0, A$ = |V>| 2 

on M. d , d € {1, 2}, which is equivalent to the nonlinear Schrodinger equation 

(1.2) i$ t + Ai> + 7 (G d (\x\) * V> - 

where G d (\x\) denotes the Green's function of the Laplacian on K d . Of physical 
interest are solutions having finite energy E and particle number (or charge) N 
given by 



(1.3) E{^) = - |V#M)r<^-^ // G d (\x-y\)\^(x,t)\^(y,t)\ 2 dxdy 
and 

(1.4) N(if>)= f mx,t)\ 2 dx, 
respectively. 

In the physical and mathematical literature the Schrodinger-Newton system in 
three space dimensions has a long standing history. With 7 designating appropriate 
positive coupling constants it appeared first in 1954, then in 1976 and lastly in 1996 
for describing the quantum mechanics of a Polaron at rest by S. J. Pekar [1], of an 
electron trapped in its own hole by Ph. Choquard [2] and of self-gravitating matter 
by R. Penrose [3]. In 1977, E.Lieb [2] showed the existence of a unique spherically 
symmetric ground state in three space dimensions by solving an appropriate mini- 
mization problem. This ground state solution u^(x), uj > is a positive spherically 
symmetric strictly decreasing function. In [4], P.L. Lions proved the existence of 
infinitely many distinct spherically symmetric solutions and claimed a proof for the 
existence of anisotropic bound states in [5] . 

While Lieb's existence proof can be easily extended to dimensions d — 4 and 
d = 5, the situation has been unclear for lower dimensions due to the lack of 
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positivity of the Coulomb interaction energy term. For the one-dimensional and 
two-dimensional problem this difficulty has been overcome recently in [6] and [7] 
and the existence of a unique spherically symmetric ground state has been shown 
by solving a minimization problem. Finally, existence and uniqueness of spherically 
symmetric ground states has been shown in any dimension d < 6 in [8] by employing 
a shooting method for the associated system of ordinary differential equations. In 
addition, for d = 1, the existence of a unique antisymmetric ground states (which 
is an eigenfunction of the parity-operator) has been shown in [6] and, for d = 2, the 
existence of unique angular excitations (eigenfunction of the angular momentum 
operator) has been shown in [7]. However, in one and two dimensions, existence 
of higher bound states remained open along with an angular momentum in the 
last case and so far only numerical studies are available indicating the existence of 
excited states, see e.g. [9]. Our main result proves the existence of such solutions 
in the attractive case 7 > 0. 

1.1. The one-dimensional case. We study the existence and uniqueness of quasi 
stationary solutions of the form 



(1.5) = m , , lim <p(\x\) = 0, 

where p = {0, 1} for an odd or even function. For solutions of the form (1.5) we have 
&(t, x) — v{\x\) and if(r),v(r) satisfy the following system of ordinary differential 
equations: 

ip" = (-fv — Oj)(f 

(1.6) v" = p 2 , r> 0. 



We suppose that v (0) are finite and v' (0) = 0. In addition is subject to the initial 
condition 

0(0) e R+, <t>'(0) = ifp = 
0(0) = 0, 0'(O)eK+ ifp=l 

The latter equation of (1.6) implies that v' > 0, and therefore for solutions 
vanishing at infinity we have u> — jv (0) > 0. Changing the variables and rescaling 
u (r) = Ar-P (r/a), V (r) = B {v (r/a) - w/7) + 1 with 



O = y/u-JV (0), A = B 



9 1 9 1 

u u 



we obtain the following system of equations: 



u " + ^P U ' = ( V - l)u 
r 



(1.7) V"=u 2 r 2 P 



subject to the initial conditions 

(1.8) it(0) = it G K+, u'(0) = 0, V(0) = 0, V'(0) = 0. 
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1.2. The two-dimensional case. We search solutions of (1.1) uniformly rotating 
with an angular velocity 0. In the rotating frame the equation (1.1) becomes 

iip t + A tp - 7$^ - O Lip = 0, A $ = \ip\ 2 

where L = —i (xd y — yd x ), the orbital angular momentum operator. 

We study the existence and uniqueness of quasi stationary solutions in the ro- 
tating frame of the form 

(1.9) ip{t,x) = ip{\x\)e- wt+vme , lim <p{\x\) = 0. 

\x\— >oo 

If, in addition, 

<p(\x\)>0, 

we call this solution the ground states for a given orbital angular momentum m, 
the other one are called excited states. For solutions of the form (1.9) we have 
<&(t, a;) = v(\x\) and <p(r),v(r) satisfy the following system of ordinary differential 
equations: 



1 m 2 
tp + - <p 2"<p = ("/v — (lj — Clm))ip 



v" + -v' = ip 2 , r>0. 
r 

We suppose that v(0), <j> (0) are finite and v' (0) = <j> (0) = 0. The latter equation 
implies that v' > 0, and therefore for solutions <j> vanishing at infinity we have 
uo — Qm — "fv (0) > 0. Changing the variables and rescaling u (r) = r~ m A <p (r/a), 
V(r)=B (v (r/a) - (u - Qm) /i) + 1 with 

a = ^Lu-nm-~fv (0), A = ^J, B=\, 

o u 

we obtain the following system of equations: 

u"+ 2{m + 1) u' = (V-l)u 
r 

(1-11) V " + 1 V > = U \^ 

r 

subject to the initial conditions (1.8). 

1.3. The general case. Equations (1.7), (1-11) are particular cases of the following 
initial value problem. 

„ 2m + d - 1 , 
u H u = (V — l)u 



(1-12) t/"^- 1 



2 2m 



subject to the initial conditions 

(1.13) u(0) =u eR + , u'(0) = 0, V(0) = 0, V'(0) = 0. 
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Here, m > may be regarded as a continuous parameter. By analyzing the solutions 
of the above initial value problem we shall prove the following result about the 
existence and uniqueness of ground states and the existence of the excited states: 

Theorem 1.1. For d = 1,2, any m > and n > the system (1.12) subject to the 
initial conditions (1.13) admits a solution (u mtn ,V m>n ) such that it m> „ has exactly 
n zeros and 



Moreover, the ground states u m fi are unique and satisfies u m fl(r) > 0, u' m (r) < 
on ]0, oo[. 

To prove the main result we use an extension of the shooting method we have 
used in [8]. Shooting methods have been successfully applied to existence and 
uniqueness of solutions in boundary value problems for second order nonlinear 
differential equations [11], [12], [13], [14] . Our paper is organized as follows: In 
Section 2 we employ a shooting method to prove the existence of ground states and 
by induction the existence of all excited states (theorem 2.8). In Section 3 study 
their decay properties to prove uniqueness by analyzing the Wronskian of solutions 
(theorem 3.4). 



We begin our study with the discussion of some general properties of solutions of 
(1.12) with initial values (1.13). Standard results will guarantee local existence and 
uniqueness of solutions, their continuous dependence on the initial values as well as 
on the parameter m and their regularity. As a consequence of local existence and 
uniqueness solutions can only have simple zeros. We shall frequently apply these 
properties in the sequel as well as the following integral equations for v! and V: 



Since uq > 0, we see from integration of the equation 2.1 that V (r) is increasing, 
since d < 2, and goes to infinity. Thus for any uq there exist a unique a > 0, such 
that V (a) = 1. This is indeed a crucial point to restrict ourselves to dimensions 
< 2. 

For the initial condition uq > of the solution (u, V) for an m > 0, we consider 
the following sets: 

Definition 2.1. 

(2.2) N m ,n '■= {uq £ K + : u has exactly n zeros} 



(1.14) 



lim u m , n {r) = 0. 



2. Existence of bound states 



(2.1) 




(2.3) 




(2.4) 



J v m 
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(2.5) M m := Q M m ,n 

n=l 

Note that since u"(0) = — u / (2m + d) all solutions start strictly decreasing. 
Therefore if u has a first critical point r\ > where u > 0, then V{r\) > 1 and 
since V is strictly increasing (see (2.1)) it follows again from (2.1) that u'(r) > 
for all r > ri . Hence every solution u ^ A/^ is strictly decreasing in the maximal 
interval (0, R) where u > 0. Obviously the sets A/^„, N m n are mutually disjoint 
and we will show that these sets forms a partition of R + . First we prove that u can 
oscillate only in a finite range and then we show that u has only a finite number of 
zero in this range. 

Lemma 2.2. Let a be such that V (a) = 1, then u has at most one zero in the 
range r > a. 

Proof. If u has two zeros for r > a, there exist a critical point r such that u (r) ^ 0, 
u (r) u" (r) < 0. Moreover from equation (2.1) we have 

u (r) u" (r) = {V (r) - l)u {rf > 

which is the desired contradiction. □ 

As consequence of the lemma 2.2, the number of zeros of u can only increase by 
one in a neighborhood of u . 

Lemma 2.3. For any u € N m ,n, there exist e > such that B(u ,e) G N m ,n U 
N m ,n+i ! where B(u ,e) is the open ball centered in u with radius e. 

Proof. We choose R > max (a, r n ) where r„ is the last zero of u. Since u has simple 
zeros it follows from the continuous dependence on initial values that for every uq 
in a open ball B (u , e), u has n zeros in the interval [0, R] and V (R) > 1. Finally 
from lemma 2.2 u cannot have more than n + 1 zeros. □ 

Lemma 2.4. u has a finite number of zeros. 

Proof. It is sufficient to prove this assertion in the range r < a. We compare 
equation (1.12) to 

„ , . 2m + d — 1 , . . 

u (r) + v! (r) +u{r)= 0, 

r 

which admits the solution: 

u (r) = w r f m + - J - J(2m+d-2)/2 (r) , 

where J(2m+d-2)/2 is the Bessel function of first kind of order (2m + d — 2) /2. By 
the Sturm's comparison theorem, see e.g. [14], u oscillate less than u in the range 
r < a. Since J(2m+d-2)/2 has a finite number of zeros in this range, the same 
properties also hold for u. □ 

The reason of the notation A/Jjf „ become from the fact that the solutions u in 
these sets goes to infinity as we will see in the next lemma. 

Lemma 2.5. For every uq <E there exist r c € (0, oo] such that lim r ^ rc \u (r) | = 

oo. 
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Proof. We have seen in the lemma 2.2 that, for large r, solutions of 1.12 have no 
zeros and, more precisely, are not oscillatory. Then, supposing lim infV^oo u (r) > 
or lim sup r ^ OCJ u (r) < leads to a contradiction. In the case liminfV^oo u (r) = 
I > 0, choosing r > a sufficiently large and integrating (2.1) gives: 

«'W - J_ (^( S ) — ^)u(s)s 2m+d ~ 1 ds + ^-J „'(f) 

> (r-f)+^V ro+1 ^W Vr>r. 

2(2m + d) v ; \r J w 

The case limsup^^ u(r) < I can be proved in the same manner. □ 

Remark 1. Since as u has only simple zeros and the preceeding lemma, we deduce 
from the continuous dependence on initial values that Af^ n , N m are open sets. 

Our main result theorem 1.1 states that each Af^ n contain at least on element. 
Obviously, the set of the ground states AT® n is nonempty if both Af m and 
are nonempty. As we will show in the proof of the main-theorem, the bound states 
are the infima of the Af m ,n sets. Thus, a necessary condition for the existence of 
bound states for each n is these infima are ever strictly positive. We will show these 
properties in the two following lemmas. 

Lemma 2.6. The set Af^o * s non em Pty- 

Proof. We want to show that u € -^mo f° r u o sufficiently large. Suppose on the 
contrary that uo € A/' m UA/'^ l for all uo > and denote ]0, Rq[ the maximal interval 
where u > and u' < 0. Let r = ^/2(2m + d). We consider the function 

/ (r) := u (r) - u 1 - — 
V r o 

It satisfied the differential equation 
9m -4-/7—1 

(2.6) /" (r) + /' (r) = V (r) u (r) + w - u (r) . 

r 

Since /' (0) = / (0) = and the right hand-side of (2.6) we conclude that / (r) > 
on [0, Rq[. Hence 

(2.7) «(r)>»o(l-|jJ on [0,i? [ 

and r < i? -Using that u is decreasing and from the equation (2.1) for V we obtain 
the bound 

V '(r) > Jili^m+i on [ 0)J R [. 

Integrating this inequality and using again that u' < yields the following estimate 
for V : 

(2.8) V ^ ^ rfiml l) r2im+1) ° n t '^ 
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We want to show that u' (ro) > provided uq sufBciently large which yields the 
desired contradiction. Using (2.7), (2.8) in equation (2.1) leads to 

i r-ra 

U 'M = -*s+3=T / {V{r)-l)u{r)r 2m+d ^dr 
r o Jo 



> 



I r u ( r ) 3 r Am+d+1 dr - — 

rl^{m + l)h U[r) 



> 



'o 



2m+d+1 (m + 1) Jo V" ^ 



T r o (m+1) jf i 1 - * 2 ) 3 t 4m+d+1 dt - 2) . 



ro \m + 

We conclude that u' (r ) > for u sufficiently large which contradicts the assump- 
tion that 7V m U7V° n = ]0, oo[. □ 

Lemma 2.7. For any N > 0, there exist uo > such that for each u < uq, the 
corresponding solution has at least N zeros i.e. u G U n >NNm,n 

Proof. First we consider equation (1.12) for uq > in the range r G [0, b], where b 
is the value at which V (b) = 1/2. We introduce a Lyapunov function F (r) defined 

by 

F(r) =u 2 {l-V)+u' 2 . 

Then F (0) = u and 

, 2(2m + d-l) ,,.2 T7 // n / n2 . n 
F (r) = -u (r) — 1/ (r) u (r) < 0. 

Hence F (0) > F (r) for all r G [0, 6] and, using V (r) < 1 we therefore have 

2u 2 > u (r) 2 . 
Then from the equation 2.1 we find 

2 

hence 

' (to + 1) (2m + d) \ 



= V(b)= / HifL? ds( ft < 6 2( m+ i) 

Jo Jo t ~ (m + l)(2m + d) 



6> 



2 U 2 



We can now compare the main equation (1.12) to the following 

„ , . 2m + d — 1 . , 1 , . 
u (r) H u(r) + -u(r) = 

V A 

which admit the solution: 

(2v / 2~\ ^ 2m+d ~ 2 ^ 2 

where J(2m+d-2)/2 is the Bessel function of first kind. By Sturm's comparison 
theorem, [14], u oscillates faster than u and has at least the same number of zero 



than J( 2m +d-2)/2 {r/V2) in the range < r < ( ( m + 1 )(| m + d ) j 2(m+1) . Therefore we 
can choose uo sufficiently small such the number of zero of u is greater or equal 
than TV. □ 
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Hence we have proved by the preceding lemmas the existence of ground states 
(see also [7]). We are now in position to prove our main existence results. 

Theorem 2.8. For any to > and neN, a m ,n hif N m ,n € ■N'mn- ^ n addition 
for every to > 0, the solution u mS) G Af^ satisfies u m ^ (r) > 0, u' m (r) < on 
]0,oo[. 

Proof. We prove this first statement by induction with the following hypothesis: 
For any n E N, a n := mij\f m _ n e Nmn- We have seen in the lemma 2.7 that 
a > 0. Since Af^.o an d ar e open sets, ao € A/^ i0 and this result is true for 
n = 0. By hypothesis a„ € A/^ „ and the lemma 2.3 applied at the point a n states 
there exist e > such that B (a n , e) c M m ,n U A/" m ,„+i. Since a„ is the infimum of 
A/" TO ,n we have the following result 

Afm,n+1 ^ and a n+l < OL n 

and it follows from the lemma 2.7 that a n+ i > 0. 

Now we suppose on the contrary that a n +\ £ Af m ,n+i- Since a„+i < a n , 
lemma 2.3 implies there exist an open ball around a n +i which is not in Af n +i in 
contradiction with the definition of a n+ \. Finally, since N^ n is an open set, it 
follows that a n +i € A/^ „• 

It remains to prove the last part of the theorem. We have seen that every solution 
u £ Af^o is strictly decreasing in the maximal interval (0, R) where u > 0. Since 
the are only simple zeros, we have in particular for the ground states: u(r) > 0, 
u'(r) < on ]0, oo[. □ 

Remark 2.9. We also prove with the preceding theorem that for all to > the 
sequence {a n }„ eN of the infima is strictly decreasing as illustrate in the following 
figure. 



The curves represente the sets Af m ,n for n = 0, 1, 2, 3, 4, 5 as a function of m in dimension d = 2. 
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3. Uniqueness of ground states 

In this section we prove that Af^ has exactly one element. First of all, we show 
that if A/"°j o had more than one element the corresponding solutions cannot cross. 
We restate the no-crossing properties of [8] in the following lemma. 

Lemma 3.1. Let 112(0) > tti(O) > and suppose that U2(r),Ui(r) exist on [Q,R] 
such that u\(r) > on [0, R\. Then U2(r) > u\(r) for all r £ [0, R\. 

Proof. We consider the Wronskian of u\, U2 defined by 

(3.1) w(r) — u' 2 (r)ui(r) — u' 1 (r)M 2 (r). 
Then w satisfies the differential equation 

(3.2) w' + 2m + d - 1 w = (V 2 - Vi)mtt 2 - 

r 

Suppose there is f £ [0, R] such that U2(r) > u\(r) on [0, f [ and u\(f) — U2(f) > 
0. Then 

w(f) = (u' 2 (f) — u[(f))ui(f) < 0. 
On the other hand we have 

V 2 '(r) V{(r) = f(A(s) ul(s))s^+ d ^ ds > 

' Jo 

on]0, f] and therefore V2 (r) > V\(r) on]0,f]. We conclude then from the differential 
equation (3.2) for w that w (r) r 2m + d — 1 j s strictly increasing on ]0,f] and since 
w(0) = we must have w(f) > which is the desired contradiction. □ 

Remark 3.2. From the no-crossing property stated in lemma 3.1 it follows immedi- 
ately that Af m , ■N'mO are intervals. More precisely, M m =]0,a[,J\f^ =]b, oo[ with 
0<a<6<oo[. Uniqueness of ground states is then equivalent to a = b. 

The important conclusion from lemma 3.1 is that two different ground state solu- 
tions cannot intersect. From the differential equation (3.2) for their Wronskian w(r) 
we see that w(r)r 2m+d ~ 1 is a nonnegative strictly increasing function. However, 
we shall prove in the sequel that w ( r ) r 2m + (i - 1 vanishes at infinity which yields the 
desired contradiction. Therefore we have to analyze the decay properties of ground 
states at infinity. 

Lemma 3.3. Let u e Af^.o- Then 

lim --V-i=l. 

1 >oo u 

Consequently, for any k £ (0, 1), 

limsupu(r)e K K vl/2(s)ds < 00. 

r >oo 

Proof. We consider the function z := —^V~^ which is well defined for all r > 
and satisfies the differential equation 



'-(•*-i)W-.(*!±^i + £) 



We also consider y := 2(m + l)U — rV which satisfies the following differential 
equation 

y' = (2m + d)V -u 2 r 2m+1 . 
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Since u is decreasing, we have from the equation (2.1) 



^M = i /' u(sfs 2m + d -'ds> 

r Jo 



u (r) 



2 



Jm+1 



2m + d 

Hence y is increasing and we have the upper bound 

2(m + l) 
V ~ r 

Now choose f such that ^n±^V' 1/2 < 1/2 for all r > f. Consider the direction 
field in the (r, z) plane for the preceding differential equation. In the set r > f, 
z > 2 we have 



z' 



> (z 2 - 1) V 1 / 2 - z V — + V- 1 ' 2 

> l - z 2 V^ 2 + l -(z + l)(z-i)V 1 ' 2 

> l - z 2 v l ' 2 (r) 
~ 2 w 

It follows that, should z(r) ever enter this region, it would blow up at finite time 
after f which is impossible. Hence z remains bounded. This also implies 

lim u\r)V- 1/2 = 0. 

r >oo 

Therefore we may apply l'Hospital's rule. We obtain 

2 ,. fu" 1 1 V \ .. (V-l 2m + d-l T __ 1/2 \ 1 

■i too r »oo\ U V 2V 6 i ) r >oo V V T ) 

This proves the first part of the lemma. 

Then for any kg (0, 1) and r sufficiently large, —^V^ 1 / 2 > k and the proof is 
completed by integrating this inequality. □ 

Now we are in position to prove our uniqueness result: 
Theorem 3.4. The set Af^o ^ as exactly one element. 

Proof. Let ui(O), 1*2(0) G a such that ^2(0) > u\(0). By lemma 3.1 the corre- 
sponding solutions wi,it 2 cannot intersect and we have U2(r) > u\(r) > for all 
r > 0. From the differential equation (3.2) for their Wronskian w(r) we see that 
ty(r)r 2m+ is a nonnegative strictly increasing function since 

(10 (r) r 2m+d - 1 ^' = (V 2 - Vi)u 1 u 2 r 2m+d - 1 > 
and w(0) = 0. On the other hand, we claim that 

lim to (r) r 2m+<,_1 =0. 

r >oo 

Indeed, since 112 cannot intersect U\ we have lim Yi{r) > V\ (r) > 1. Trivially, 

r— »oo 

i)(2m+d) ■ From the integral equation (2.1) for u^i 

u' 2 (r) r 2m+d - 1 = /""(^ 2 («) - l)u 2 {s)s 2m+d ~ 1 ds 
Jo 
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and the decay properties of u 2 given in lemma 3.3 it follows then that u' 2 r 2m+d 1 
and U2r 2m+d ~ 1 are uniformly bounded. Therefore 

\w{r)r 2m+d - 1 \ < \u l \\u' 2 r 2m+d - 1 \ + |u / 1 ||u 2 r 2,n+<, - 1 | < ci|ui| + c 2 \u[\ 

for some positive constants Ci , c 2 which concludes the proof. □ 
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